Variable size particle probing of the DLA cluster properties 
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We develop a technique for probing harmonic measure of the diffusion limited aggregation (DLA) 
cluster surface with the variable size particle and generate one thousand clusters with 50 million par- 
ticles using original off-lattice killing-free algorithm. Taking, in sequence, the limit of the vanishing 
size of the probing particles and then sending the growing cluster size to infinity, we achieve the un- 
precedented accuracy determining the fractal dimension, D = 1.7100(2) crucial to characterization 
of geometric properties of the DLA clusters. 
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Kinetic interfaces that evolve into two dimensional 
fractals via the diverse stochastic growth processes are 
ubiquitous to nature. The systems and physical pro- 
cesses that exhibit kinetic roughening and fractal struc- 
ture range from firefronts and bacterial colonies to den- 
drites of various nature, to domain walls in magnets 
and ferroelectrics, to liquids penetrating porous media, 
and many others. The study of two-dimensional frac- 
tals associated with the kinetic roughening is an exciting 
and mature branch of statistical physics - see excellent 
works P, 0] for exhaustive reviews. The generic process 
governing a good part of these phenomena is the so called 
two dimensional aggregate growth, which is commonly 
modeled as the diffusion limited aggregation (DLA) Q 
and its generalization DBM 0| well capturing most of 
the properties of the random aggregates 0- Past ana- 
lytical and numerical studies advanced impressively our 
understanding of DLA (see again 0, 0), yet there are 
still several critical issues that remain unresolved, like the 
controversies related to the multiscale and fractal nature 
of DLA, to name a few. One of the central questions is 
the precise value of the DLA fractal dimension D, the 
knowledge of which is crucial for full description of geo- 
metric properties and characterization of the DLA clus- 
ters. In particular, the exact value of D is necessary for 
evaluating DLA lacunarity at large growth times. 

In the planar, d=2, geometry the analytical results (see 
for example [3) suggest that the DLA fractal dimension 
lies in the range 15=1.67—1.72; certain rational values 
like D=5/3, in or D^n/IQ, in were also pre- 
dicted. On the other hand, Mandelbrot Q argued that 
DLA clusters fill the whole space, i.e. that D=2. The 
direct simulations usually produce 13=1.715(4) for the 
DLA fractal dimension I9I I1OI: the simulations |ll| based 
on the conformal mapping 

0,13 yield D^l. 713(5). 

Conventionally the fractal dimension D is determined 
from the fit of the growing aggregate (cluster) size, R, to 
the functional dependence 
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where N is the number of particles in the cluster. The 
cluster size, R, can be defined as the radius of the de- 
position Rdep—{ri), where is the position of the i-th 
particle, and brackets stand for the averaging over the 
ensemble of clusters. The accuracy achieved in the deter- 
mination of D in the pas t publications is characterized by 
the fluctuations J^d= {{D^)-{D)^) /{D)^- The 

latter decays as !Fd(xN~°-^^o:^N^^^^ (rather than naively 
expected ocl/iVlniV). As a matter of practice, this im- 
plies that in order to achieve the next level of accuracy 
[one more (forth) digit] one is required to consider clus- 
ters with the number of particles exceeding by the factor 
of 1000 those used in simulations of Refs. |9l ll0lllll |. The 
clusters of such a size, about A^=10^, are not accessible 
via standard approaches due to the computational speed 
and computer memory limitations. 

In this Letter we develop and report on the alternative 
approach allowing to reach higher accuracy of the frac- 
tal dimension measurement. We exploit the fact that in 
simulations the diffusing particles always have some finite 
physical size 6 and that while the true fractal dimension 
D should not depend on S, the results of simulations 
appear to depend on the size of the particles probing 
harmonic measure. Our idea is to use this dependence 
in order to parameterize the fractal dimensionality as a 
function of S and, upon finding D{6) in a series of nu- 
merical experiments, define the fractal dimension as 
D= lim5_,o D{S)- We will show that our approach allows 
for a dramatic improvement in the precision in determin- 
ing D. 

The accessibility of a particular set of the points at the 
cluster interface is characterized by the probabilities pk 
for a particle diffusing from the infinity to hit the cluster 
at this given subset, Ffc, of the surface. For the interface 
of the form of an ideal circle (and the symmetric diffu- 
sion) , the probabilities are distributed uniformly, and the 
probability density is just 1/2tt (the harmonic measure). 
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FIG. 1: Fractal dimension with error bars as function of the 
cluster size, computed with Rdep- Solid line connecting circles 
is guide for the eye. 

In an irregular cluster most of the probability density is 
concentrated around the cluster's tips, whereas the fjords 
are screened (this effect is often referred to as the "tip ef- 
fect" and/or the "Faraday screening" |l,[ll[Il[I3). The 
conventional approach to probing the harmonic measure 
is as follows: (i) The DLA cluster containing particles 
is generated; (ii) The positions where the next M par- 
ticles hit the border of the cluster for the first time are 
stored but the particles do not remain attached to the 
cluster; (iii) The deposition radius is calculated as 

1 *^ 

i?<iep(iV)«-5^n, (2) 

i=l 

and the sum converges to the integral over harmonic mea- 
sure Rdep{N)= Jdqr for sufficiently large M. 

Shown in Figure ^ is the typical variation of D calcu- 
lated using relation Q with Rdep, and representing the 
half period of oscillations of fractal dimension D with 
the size of the cluster. The oscillations are related to 
the spatially nonuniform growth of the cluster: different 
branches grow with the different speed, and the new- 
born sub-branches often outgrow the parent branches. 
The amplitude of the oscillations decrease slowly with 
the cluster size thus complicating precise determination 
of the value of D. 

Figure El shows the part of the cluster branch. The 
spots where the probe particles hit the surface are marked 
by the color with the intensity proportional to the log- 
arithm of the particular probability pk for the particle 
to hit the segment Ffc. For the sake of better visualiza- 
tion we choose the length of segments equal to the one 
pixel of the figure. Uncoupled colored segments, Ffc (see 
Fig. 121) , constitute only the part of the surface accessible 
to the diffusing particles. The total length of the accessi- 
ble cluster interface depends on the size of the probe par- 




FIG. 2: (Color online) Fragment of the cluster with iV=10* 
particles and positions of hits by probe particles with radius 
S—1. Color is proportional to logarithm of the hit probabil- 
ities, denoted by the color bar at the top. Probabilities are 
larger near the tips and smaller inside the fjords. 

tides. We choose the probe particle size (5 as a parameter 
of the particular measurement. The value of S controls 
two processes: (i) the 'geometrical' process which is the 
penetration of the particles inside the fjords where the 
geometrical bottleneck of the fjords may prevent [TtIIi^ 
sizable particles to go through, and (ii) the probabilistic 
process of screening; the latter becomes more effective 
with the growth of the particle size, S. 

It is important to stress the difference between the 
measurement of the length of the fractal surface jl^] and 
the measurement of the harmonic measure on the sur- 
face within our approach. In the former case, the total 
length of surface grows with the decreasing scale division 
value 10 (i.e. the length of the subsets Ffc are equal to 
the ruler size). In the latter case, the probability pk for 
the particle to hit the subset F^, saturates as S^Q. 

Shown in FigureOis the number of the surface particles 
having been touched by the probe particles as function of 
the probing particle size. The solid line is given by the ex- 
pression Nreach—Nsurf/il+S/So)", with Q;=0.69(2) and 
So—2.2{l). The limit of vanishing size S^Q gives the lim- 
iting number of the reached particles as the total number 
of surface particles, 7Vs„r/~14, 662, 525 (out of total 20 
millions in the cluster). 

The dependence of D upon N for different 5 is shown 
in Fig.^ As we have already mentioned, the dependence 
on A^ within the given interval of A^ and for the fixed 6 
is not monotonic for <5<3. Thus, we fit results by the 
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FIG. 3: Number of reachable sites of the cluster as function 
of the probe particle radius (circles) with the fit described in 
the text (solid line). 
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FIG. 4: Fractal dimension (symbols) as function of A'^ for 
S =0.1, 0.3, 1, 3, 10, 30, 100 from bottom to top. Solid lines 
connecting circles is guide for the eye. 



D{S; N) = D(N) + AS^^ 
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and take the limit of 5^{) for the fixed values of the clus- 
ter size N . The resulting values of D{N) are shown in 
Fig.|Sl Note that now the dependence D{N) has become 
monotonic. This can be understood as the result of the 
effective averaging the competition between the growing 
branches out. We have examined the ensemble of 1000 
clusters with 50 million particles for each value of S. The 
number of particles M in each event of the measurement 
necessary to achieve the desired accuracy in Rdep of about 
0.1%, was typically several tens of thousands. We have 
found no difference in results when using larger number 
of probe particles. The error bars represent combined 
errors from both the ensemble average and the fit. Thus 
from Fig. [S] one finds at the end of the day the ulti- 
mate value of the fractal dimension, D=l. 7100(2). This 
is an unprecedented accuracy in the measurements of D 
exceeding by the order of magnitude the results known 
from the literature. 

It is instructive to use the same ensemble of clusters in 
order to evaluate the fractal dimension via the traditional 
approach introducing the deposition radius as Rdep— 
the mean square displacement i?2 = \/ (r^), and the radius 

of gyration Rg = \J 'J2iLi{'''^) ■ The results of the fit to 
the form where R stands for one of the above radii, 
are presented in Table |I| Note, that all the obtained 
values of D are larger than that of L)=l. 7100(2) derived 
by the method of present work. 

Since DLA clusters grow randomly, a center of the clus- 
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the plane. Accordingly, the distance from the original 
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FIG. 5: Fractal dimension as function of A'^ for the limiting 
value of (5 = 0. It reach value D = 1.7100(2) (plotted with 
the dotted line) in the limit of large cluster size A*'. 



position of a seed particle to the center of mass grows as 
\/N and the number of particles is proportional to the 
time of the random walk ji^l • The average value of Rm 
for iV=5T0^ is (i?M)~1000 which should be compared 
to (i?g)«30000 and the penetration depth ^«7000. The 
average angular position of the center of mass (averaged 
over an ensemble) is also a stochastic function of time, 
although at each realization of the cluster angular cor- 
relations are observed (some of the branches grow faster 
within a given time interval |20j'). 

The important question now is whether the choice of 
the coordinate frame influences the final result and the 
value of the fractal dimensionality. Indeed, when deter- 
mining D one can choose the origin at either (i) the po- 
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r to seed 


r to center-of-mass 


Rdep 

R2 

Rg 


1.7098(12) 

1.71155(56) 

1.71149(30) 


1.7111(6) 

1.71149(54) 

1.71133(30) 



TABLE L Fractal dimension estimated by the standard 
method using the reference frames centered at the seed parti- 
cle (left column) and at the center of gravity (right column) 
respectively. 



sition of the seed particle, or (ii) at the (evolving with 
time) position of the center of gravity, or else (iii) at the 
'center of the charge gravity' (the latter is most appro- 
priate for the DMB case). We have performed averaging 
according to Eq. |2Jl finding D from the fit to Rdep, R2, 
and Rg, placing the origin of the reference frame to the 
seed particle (left column of the Tabled and to the cen- 
ter of gravity (right column). One sees that for all the 
quantities the choice of the reference frame is irrelevant 
(within the accuracy of the computation), although the 
ensemble of clusters should be large enough to insure this 
convergence. 

In conclusion, we have developed a technique for the 
high precision analysis of the geometrical properties of 
DLA clusters, in particular the evaluation of its lacunar- 
ity in the long time limit. Our approach offers a perfect 
tool for further advance in our understanding of the fjord 
screening behavior. In particular, the long standing prob- 
lems of the behavior of the Pmim which is the minimal 
growth probability (in DBM language) or probability to 
hit given segment of the cluster surface (in DLA lan- 
guage) can be efficiently addressed via the developed ap- 
proach. The behavior of the latter probability is related 
to the phase transition in the multifractal spectrum \2l\ . 
To the best of our knowledge, this probability was esti- 
mated only within the "tunnel configuration" |l6| tech- 
nique and had never been measured in the simulations of 
the "typical configuration" 22]. 

We are pleased to thank S. Korshunov for important 
and enlightening discussions. This work is supported by 
the US DOE Office of Science under contract No. W31- 
109-ENG-38 and by the Russian Foundation for Basic 
Research. 



[2] A. Bunde and S. Havlin, eds.. Fractals and Disordered 

Systems (Springer, Berlin, 1996). 
[3] T.A. Witten and L.M. Sander, Phys. Rev. Lett. 47, 1400 

(1981). 

[4] L. Niemeyer, L. Pietronero, H.J. Wiesmann, Phys. Rev. 

Lett., 52, 1033 (1984). 
[5] H. Gould, F. Family, and H.E. Stanley, Phys. Rev. Lett. 

50, 686 (1983). 
[6] M. Muthukumar, Phys. Rev. Lett. 50, 839 (1983). 
[7] M.B. Hastings, Phys. Rev. E 55, 135 (1997). 
[8] B. Mandelbrot, Physica A 191, 95 (1992). 
[9] S. Tolman and P. Meakin, Physica A 158, 801 (1989); 

Phys. Rev. A 40, 428 (1989). 
[10] P. Ossadnik, Physica A 176 (1991) 454, ibid. 195, 319 

(1993). 

[11] B. Davidovitch, A. Levermann, and I. Procaccia, Phys. 

Rev. E 62, R5919 (2000). 
[12] M.B. Hastings and L.S. Levitov, Physica D 116, 244 

(1998) . 

[13] B. Davidovitch, H.G.E. Hentschel, Z. Olami, I. Procac- 
cia, L.M. Sander, and E. Somfai, Phys. Rev. E 59, 1368 

(1999) . 

[14] T.A. Rostunov, L.N. Shchur, JETP 95, 145 (2002) 
[15] A. Yu. Menshutin, L. N. Shchur, Phys. Rev. E 73, 011407 
(2006). 

[16] R. Blumenfeld, A. Aharony, Phys. Rev. Lett. 62, 2977 
(1989) 

[17] R.C. Ball, R. Blumenfeld, Phys. Rev. A 44, R828 (1991). 
[18] H.G.E. Hentschel, Phys. Rev. A 46, R7379 (1992). 
[19] B.B. Mandelbrot, Science 155, 636 (1967). 
[20] A.Yu. Menshutin, L.N. Shchur, and V.M. Vinokur, un- 
published. 

[21] J. Lee, H.E. Stanley, Phys. Rev. Lett. 61, 2945 (1988) 

[22] J. Lee, P. Alstrom, and H.E. Stanley, Phys. Rev. Lett. 
62, 3013 (1989) 

[23] Our definition of the fractal dimension does not depend 
upon the size of the particles the cluster itself is composed 
of. Therefore, we let the size of the particles in the cluster 
to be unity and measure the size of the probing particles 
in the units of the cluster particles. 

[24] There are some attempts in literature of extracting the 
fractal dimension of the cluster using fit RmocN^^'^ . 
However no justification for such fit is available to the 
best of our knowledge. One rather expects the center of 
cluster mass to perform a conventional random walk. 



[1] T. Halpin-Healy and Y.C. Zhang, Physics Reports 254, 
215 (1995) 



